S-wave contact interaction problem: A simple description 



Shina Tan 

James Franck Institute and Department of Physics, University of Chicago, Chicago, Illinois 60637 

The s-wave contact interaction problem has a very simple structure and a simple and straightfor- 
ward description. This is not a standard paper, and technicalities are avoided. A bare and simple 
picture of the problem is presented, in the hope that it will be helpful to our theoretical study of 
the problem. 
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The s-wave contact interaction, or s-wave resonant in- 
teraction, is purely a quantum mechanical effect, and has 
no classical counterpart, unlike the Coulomb interaction. 
This concept (s-wave contact interaction) however has 
such a broad range of applications that numerous theo- 
retical and experimental papers are published, which are 
closely related to it. In this regard, it is not unlike the 
Coulomb interaction. 

Two nonrelativistic particles in three-dimensional 
space are said to have s-wave contact interaction if they 
satisfy the following conditions: 

• they have no direct interaction, if their distance or 
relative orbital angular momentum is nonzero; 

• the s-wave component (ie, spherically symmetric 
component) of their relative wave function obeys a 
simple linear boundary condition at zero distance, 
as below. 

For simplicity, we consider particles of equal mass m, 
and use the natural unit h = 1. Generalizations to un- 
equal masses and/or restoration of H is easy. 

We first omit the center-of-mass degree of freedom, 
and concentrate on the relative motion. The Schrodinger 
equation is 
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where ip = ip(r, t), r is the relative wave vector of the two 
particles, and t is time, r = |r| is the distance between 
the particles. 

At short distances r — > 0, we can omit the energy, and 
obtain an even simpler equation, — V 2 t/> = 0. Different 
partial waves, associated with different angular momen- 
tum quantum numbers I, are decoupled: ip = X^o^' 
such that 
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and thus ipi — Ai(v)r + Bi(r)r~ , where f = r/r is the 
direction of r. Since there is no interaction in channels 
of nonzero angular momenta, -Bz(r) = for I > 0. In the 
s-wave channel, however, the interaction is denned by a 
linear constraint on Aq and Bq: 
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and here enters the scattering length a; we assume it is 
a nonzero constant; it is either positive or negative, or 
even infinity (the so-called unitarity limit), a is the only 
parameter characterizing the s-wave contact interaction. 
Equation J2J), as simple as it is, is of paramount impor- 
tance in the s-wave contact interaction problem. Note 
also that the s-wave contact interaction can not exist be- 
tween identical fermions; for other particle pairs (includ- 
ing a pair of fermions at different internal states), such 
interaction is allowed. 

Now the wave function is generally of the form 



tf>(T,t) = B (t) [--- 
r a 



+ 0{r) (r^0). (3) 



To facilitate the study of much more complicated cases, 
in which a few or many particles participate in the s-wave 
contact interactions, we introduce the short-range selec- 
tors A(r) and l(r) [l| [this l(r) happens to involve the 
same English letter as the angular momentum, unfortu- 
nately], whose role is to selectively pick up the short- 
range behavior of the wave function (or some other func- 
tions that may arise in the problem) . These selectors are 
possible if we have the general notion that mathemati- 
cal tools (as well as any physical concepts) can be freely 
reshaped to suit our needs: 

• truth value (namely the proposed theory, in its 
range of validity, must not contradict our experi- 
ence); 

• coherence and simplicity (Ockham's razor); 

• ease to handle (namely the logical path from the- 
ory to experience should be shortened as much as 
possible). 



Here are the basic properties of these selectors: 
A(r) = l(r) =0 (r ± 0), 
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J d 3 rA(r) =1, fd s r£± =0, (4b) 
d 3 r/(r) =0, /d 3 r| =1, (4c) 



J d 3 rA(r)f = J d 3 r/(r)f = 0. 



(4d) 
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These selectors have momentum representations: 

A(k) =/d 3 rA(r)exp(-ik-r), 
L(k) = Jd 3 W(r)exp(-ik-r), 

and they have properties 



A(k) = 1 (|k| < oo), 
d 3 fc A(k) 
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A(-k) = A(k); 
L(k) = (|k|<oo), 
d 3 fc L(k) _ 
(^"fc 2 " ~ ' 
L(-k) = L(k). 

These equations are not contradictory, because 
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in contrast with the traditional notion about integrals. 

These two selectors span a linear space, called the se- 
lector space. 

Finally, we define the important ^-selector, which lies 
in the selector space, 



r}(r) = A(r) + 
7?(k) = A(k) -f 
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and it satisfies the important equation 

= 0, 
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but for functions that are nonsingular at r = 0, rj(r) 
behaves like the delta function. 

Equation © can now be written as 



d 3 r?7(r)V>(r,i) = 0, 



(9) 



which will play a crucial role in the s-wave contact inter- 
action problem. 

We now return to Eq. (JIJ. Because the wave function 
has a singular term Bq jr, the kinetic energy operator pro- 
duces a term AttBqS^ (r)/m, which has to be canceled, 
if we want to make the Schrodinger equation valid for all 
r. The simplest way to do this is to use the L-selector, 

*JU(r,*) = -—#",*) " -S {3) (r) f dVZ(r')V(r',i), 
at to m J 

(10) 



which must be solved in conjunction with Eq. @ in order 
to determine the wave function. 

There is also a simple way to combine the above equa- 
tions 

d V 2 

+ 1,5(3) (r) f d 3 r / + C ^ (r /)] ^ t) 
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whose solutions necessarily satisfy both Eq. JSJl and 
Eq. H10() . Here C is any nonzero length. 

The selectors —L + (r), for all values of f, form a 
straight line in the selector space, and I call it Olshanii- 
Pricoupenko straight line. Simple mathematics shows 
that Eq. (1) of Ref. [2J exactly corresponds to this line. 

Note that our A-selector is a completely different thing 
from the momentum scale "A" in 0/ or in any other 
papers. 

The momentum scale "A" in is related to the length 
£: "A"= 1/o-4tt/C. 

The Olshanii-Pricoupenko pseudopotential is now ex- 
pressed in terms of the selectors as 

V*(r, r') = -5^ (r) [-Z(r') + Cv(r')] , (12) 
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and, in particular, when £ = 47ra, we get back to the 
Fermi pseudopotential, which can now be expressed as 



V(r,r') 
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whose momentum representation is simply 
V(k,k') = — A(k'). 
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It was first noted by Olshanii-Pricoupenko that the 
Fermi pseudopotenial is a special case of many equivalent 
choices. 

Such equivalence is due to the fact that the wave func- 
tion satisfies the constraint Eq. 

The new insights that we have are: 

• the short-range behavior of the wave function is 
confined in a linear space (and is subject to an ad- 
ditional linear constraint); 

• we can construct its dual linear space, called the 
selector space, whose elements are called the selec- 
tors; 

• each selector corresponds to a well-defined general- 
ized function, and transformation between different 
representations (eg, between coordinate space and 
momentum space) is easy and straightforward; 

• this allows us to formulate the whole problem in a 
simple, unified, and flexible way, and significantly 
facilitates our study of the whole system; 
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• the old forms of the pseudopotentials in terms of 
partial derivatives are both inelegant and not flex- 
ible (for example, no one was ever able to trans- 
form them to the other representations), and are 
no longer needed, once the selectors are clearly de- 
fined. 

We now proceed to a system of fermions in two "op- 
posite" "spin" states, subject to some external potential 
V^ x t; the s-wave contact interaction is present between 
the opposite spin states. The system is described by two 
equations: 



d 3 r ?7(r)^ T (R + r/2)V>j.(R - r/2)|0> = 0, (14a) 

E = J2 I d 3 r^t( r )[_v 2 /2m + K !X t(r)]^(r) 
_JL f d 3 rd 3 r> T t (r)?A{(r)?A i (r-r72)V' T (r + r72)/(r'), 
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(14b) 

where \<f>) is any physically allowed state. We can of 
course combine the two equations as previously. 

It is shown in pj that the energy of the system can 
be exactly expressed in terms of the one-particle reduced 
density matrix only: 



E = J2^ k ^ n ^ + J2 J d 3 rKxt(rK(r), 

kCT (7 



(15) 



where n^ a and p a (r) are the momentum distribution 
and number density distribution, respectively. Note that 
7/(k) = +1 for any finite k, but the peculiar behavior 
of ry(k) upon integration allows the existence of negative 
energy states, including the two-body bound states. 

Equation i|14a[) directly leads to many identities that 
must be satisfied by various reduced density matrices; see 
for details. 

The same approach can also be used to study the Bose 
gas with s-wave contact interaction. However, the same 
level of completeness will not be reached until we also 
treat the Efimov effect, inherent in this system, clearly. 
This issue was not noticed in |2|. 

Generalization of this formalism to two dimensions is 
easy. The relative wave function is of the form B$ ln(r/a) 
plus higher order terms at short distances. We can there- 
fore define two selectors rj(r) and l(r) such that 



d 2 r?y(r) =1, J d 2 r^(r)(27r)- 1 \n{a/r) =0, 
d 2 r/(r) =0, Jd 2 r/(r)(2 7 r)- 1 ln(a/r) =1, 

Their Fourier transforms 



and A(r) = Z(r) = for r ^ 0. 
can be shown to satisfy: 



|k|>K 



d 2 fc ??(k) 
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= J- m ^J-lnH^, (16b) 
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L(k) = (|k| < oo), 



d 2 fc L(k) 
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?7(k) = 1 (|k| < oo), 



(16a) 



where 7 = 0.5772 • • ■ is Euler's constant, and K is any 
positive momentum scale. 

Now the two-component Fermi gas in two dimensions 
with s-wave contact interaction satisfies two equations 
which look formally identical with Eqs. (|14afl and l|14b|l . 
except that the three-dimensional spatial integrals are 
replaced by two-dimensional ones. 

Using the logic similar to that in Q , we can show that 
the two-dimensional system satisfies the energy theorem 



£ = E^ k )|^ + E / d 2 rV cxt (r) Pa (r), (17) 

which looks almost identical with its three-dimensional 
counterpart, except that now the meaning of ?7(k) is dif- 
ferent. 

This formalism is much more straightforward than the 
one used in Q. Also the energy theorem is unknown to 
people before. 

As soon as we modify the traditional notion about in- 
tegrals, in the way demonstrated in [| and here, the 
ultraviolet divergence problems suddenly disappear. 

While this approach is certainly not the first one which 
removes ultraviolet divergence problems, a lot of advan- 
tages exist in this formalism. Let us compare it with two 
other approaches: 

• the traditional pseudopotential approach, which 
has been further developed by Olshanii and others 
□ 

• the bare coupling constants plus dimensional regu- 
larization plus renormalization approach. 

The first approach never gained the flexibility to reach a 
simple momentum representation. The second approach, 
confined by the traditional notion of integrals, involves 
complicated intermediate steps (in which the dimension- 
ality of the system has to be changed), divergent bare 
coupling constants, and subtraction procedures after the 
rcgularizations, and is thus a long detour to a simple 
problem. 

The s-wave contact interaction problems, in both three 
and two dimensions, are now formulated easily with the 
short range selectors, and other known approaches can 
hardly compete with it in this area. The author however 
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does not exclude the possibility of the appearance of a 
still simpler new approach. 

Whether the ideas here can be further extended to 
other important systems, in particular relativistic quan- 
tum field theories (in which the divergent bare coupling 
constants plus dimensional regularization plus renormal- 
ization approach currently dominates), remains an open 
question for the moment. We will get an answer to this 
question, either positive or negative. 

Finally, let us discuss a little bit about some ideas un- 
derlying the present approach. Since the prospects of this 
approach beyond the current two systems are still not 
clear, any discussion here should be regarded as rather 
tentative. 

The traditional notion about integrals can be illus- 
trated with a simple example, J d 3 k/k 2 , which contains 
divergent contribution from large k. Here k is either mo- 
mentum, or spatial vector, or even some other quantity. 
Let us assume it is momentum from this point on. This 
apparent divergence seems to stem from an assumption, 
namely that any short- distance physics, including Planck 
scale physics, and maybe even beyond Planck momen- 



tum, is really there, actively underlying that cloud of 
ultracold atoms, even when we have not fed any particles 
with very large energies. An even deeper assumption un- 
derneath this one, is reductionism: short-range physics 
completely determines long-range physics, and the for- 
mer is more fundamental than the latter. Whether these 
assumptions are valid, is of course a question beyond the 
scope of this tiny article, and we shall not attempt to 
discuss it further here. 

The intuition that we have gained from our new gen- 
eralized functions is: short-range physics emerges when 
we accelerate particles to its associated energy scale; 
whether it is really there, before we accelerate particles 
to that energy scale, seems to be a philosophical ques- 
tion - neither provable nor disprovable. However, in the 
present approach, the short-range physics appears to be 
somehow activated or maybe, more strangely, created, af- 
ter the particles are raised to high energies. Is this picture 
an illusion caused by this mathematical method, or does 
it suggest some element of truth? 

The author is indebted to Z. Q. Ma and C. Chin for 
some discussions. 
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